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A system of Na atoms of n-levels interacting dipolarly with £ modes of electromagnetic field is 
considered. The energy surface of the system is constructed from the direct product of the coherent 
states of U(n) in the totally symmetric representation for the matter times the £ coherent states 
of the electromagnetic field. A variational analysis shows that the collective region is divided into 
£ zones, inside each of which only one mode of the electromagnetic field contributes to the ground 
state. In consequence, the polychromatic phase diagram for the ground state naturally divides itself 
into monochromatic regions. For the case of 3-level atoms in the H-conhguration in the presence 
of 2 modes, the variational calculation is compared with the exact quantum solution showing that 
both are in agreement. 
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I. INTRODUCTION 

Research in the interaction between matter and elec¬ 
tromagnetic radiation has proved very fruitful no less for 
its interest in quantum information theory, a field emerg¬ 
ing from quantum optics and information science. In 
order to study the interaction between matter and radi¬ 
ation a simple but non trivial model, which considers Na 
identical atoms of 2-levels under dipolar interaction with 
one-mode of an electromagnetic field was proposed by 
Dicke [1]. The Tavis-Cummings model [2, 3] is a simplifi¬ 
cation which considers the rotating wave approximation 
(RWA), eliminating the terms in the Hamiltonian that 
do not preserve the number of excitations in the sys¬ 
tem. The particular case of a single 2-level atom with the 
RWA approximation is known as the Jaynes-Cummings 
model [4] and has analytical solution. For the case with¬ 
out the RWA approximation, Hepp and Lieb found in the 
early seventies that the system suffers a transition from 
the normal to the superradian phase [5, 6] 

Whereas n-levels atoms are only an approximation to 
real atoms, the design and construction of artificial quan¬ 
tum structures allows one to refer to the so-called arti¬ 
ficial atoms [7-9] that possess a finite number of levels. 
Hence it is interesting to consider n-level systems without 
loss of generality. 

A first generalisation of the models involving 2-level 
atoms is the consideration of atoms of three or more lev¬ 
els [10-12] which, as is well known, possess different ar¬ 
rangements: the S-, V- and A-configurations for 3-level 
atoms, and for 4-level atoms the configurations S, A, y, ■0, 
N, .... Clearly, systems of three or more levels permit 
to consider dipolar interactions with one, two or more 
modes of electromagnetic field. In fact, depending on 
the atomic configuration, there is a maximum number 
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of different modes that may produce atomic transitions 
via the dipolar interaction. For 3-level atoms, the in¬ 
teraction with one-mode of the electromagnetic field has 
been extensively studied [10, 11, 13-15]. More recently, 
treatments including the diamagnetic term show that, 
in spite of the Thomas-Reiche-Kuhn rule for multilevel 
atoms, there is a broad range of physical interaction pa¬ 
rameters for which it is possible to have a superradiant 
transition [16]; and under a degenerate parametric am¬ 
plification process and with an external laser driving field 
the system exhibits the effects of spontaneous emission 
of atoms and cavity loss in the phase diagram [17]. Two- 
modes of electromagnetic field have also been considered 
in [10, 18, 19[. In particular, the phase diagram in the 
A-configuration was shown to be constituted by a normal 
and two superradiant phases [20, 21]. 

In this work we study the properties of the variational 
ground state of Na atoms of n-levels interacting dipolarly 
with £ modes of electromagnetic field, where each mode 
promotes transitions only between two given atomic lev¬ 
els and where £ is the maximum number of modes that 
are permitted by the dipolar interaction. The variational 
ground state is evaluated by considering as test function 
the direct product of Heisenberg-Weyl coherent states for 
the field and U(n) coherent states for the matter. The 
variational energy surface is evaluated as the expectation 
value of the Hamiltonian for this test state. We find in 
general that critical points both at zero and at infinity 
(for the matter variables) provide minimum values of the 
energy, and that taking particular values of them (at zero 
or at infinity) reduces the system to (n — l)-levels, which 
may again be reduced to an (n—2)-level system and so on, 
through this iterative process in the end one has reduced 
systems of the well-known Dicke model, which in the vari¬ 
ational method can be solved analytically. The method 
is exemplified for the cases of 3- and 4-level atoms. Gen¬ 
eral expressions of the critical points that provide the 
minimum energy surface are obtained as functions of the 
dipolar intensities, and from them the corresponding nor- 
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mal and collective (superradiant) regions are determined 
for the case of n-level atoms. 

We show that the collective region is divided into I 
different zones where in each zone only one-mode of elec¬ 
tromagnetic field contributes to the ground state, i.e., the 
polychromatic phase diagram is divided into monochro¬ 
matic regions, and crossing from one to the other rep¬ 
resents a first order phase transition. These transitions 
manifest themselves as discontinuities in the atomic pop¬ 
ulations. On the other hand, crossings from the normal 
region to any of the collective monochromatic region give 
place to first or second order phase transitions; those of 
second order correspond to bifurcations, while those of 
first order form Maxwell sets [22]. This analysis leads 
to a universal relationship between expectation value of 
the number of photons in the mode and the 

quadratic fluctuations of the number of atoms in the level 
j, {AAjj)'^. We propose this relation as an experimental 
criterion to detect the transition between the normal and 
superradiant regimes. 

This paper is organized as follows: Section II describes 
the Hamiltonian for a system of Na atoms of n levels un¬ 
der dipolar interaction with £ modes of electromagnetic 
field, where only one mode promotes transitions between 
two given atomic levels. In section III, the variational 
ground state (test function) is defined for the field and 
matter contributions. The energy surface is calculated 
and the method (which involves an iterative reduction of 
the system) to find the critical points that minimize the 
energy is discussed. In section IV the method is exem¬ 
plified for the case of 3-level atoms interacting dipolarly 
with 2-modes of electromagnetic field is considered for 
the three atomic configurations. The minimum energy 
surface, the order of transitions, and the statistical prop¬ 
erties of the variational ground state are established. In 
section V the critical points of the energy surface for the 
case of 4-level atoms in dipolar interaction with i modes 
{£ = 3,4) are obtained by means of the iterative proce¬ 
dure. Section VI generalizes the expressions to the case 
of n-level atoms interacting dipolarly with i modes. Also, 
we discuss in general the quantum phase diagram, the or¬ 
der of the transitions, the expectation values of field and 
matter observables and the corresponding results when 
the RWA approximation is considered. Finally in sec¬ 
tion VII the numerical exact quantum solution is given 
for 3-level atoms in the S-configuration for Na = 1 and 
Na = 2 particles. In section VIII we give some conclud¬ 
ing remarks. 


II. MODEL 

We assume that any given pair of atomic levels ojj and 
Wfc can be connected by at most one radiation mode de¬ 
noted by the photon annhilation ajk and creation 
operators. The Hamiltonian describing a system of Na 
atoms of n-levels on dipolar interaction with t modes of 


electromagnetic field is then written as [15] 

H = HD + H,nu ( 1 ) 

with a diagonal contribution 

n n 

He, = ajk + tUj Ajj . (2) 

j<k j=l 

We have taken Ti = 1, and the collective atomic opera¬ 
tors Akj obey a U(n) algebra satisfying the commutation 
relations 

[A^Tn; A^j] = Sjnk -£^lj ^jl -^km ; (3) 

with the first order Casimir operator determining the to¬ 
tal number of atoms Na given by 

n 

Y,Auk = Na, (4) 

fc=i 

which is conserved. We denote field frequencies by ££ljk 
and assume that the atomic frequencies satisfy uJi < uj 2 < 

* ■ ■ ^ UJji- 

The interaction Hint between the atoms and radiation 
field involves contributions of the form {Ajk-£Akj){a,jkA 
a^j.) representing the fact that a transition between the 
levels j and k is only promoted by the mode Thus 
we may write 

1 " 

Hint = -J= £Ljk {Ajk + Akj) yljk + ®j/c j (5) 

where fijk denote the dipolar intensities. 

The maximum number of dipolar interaction strengths 
of an n-level system is imax = n{n — l)/2 — (n — 2); of 
course the real number £ of modes present in the system 
will be less than or equal to £max and it depends of the 
considered atomic configuration with its allowed dipolar 
transitions. 


III. ENERGY SURFACE 

In order to find the minimum energy surface, we pro¬ 
pose as a variational test function a coherent state con¬ 
taining contributions of both field and matter. 

The coherent field contribution is well known, 
since each electromagnetic mode £ljk is described by 
the Heisenberg-Weyl state \ajk) obeying ajk \c(jk) = 
ajk Ictjk)- The dimension of vector a = {ajk), j < k = 
1, 2, ..., n depends on the number of allowed dipolar 
transitions in the given atomic configuration. The total 
coherent electromagnetic field contribution is then given 
by the direct product of these states, i.e., 

71 

|a) = 0 \^jk) ■ (6) 

j<k 
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The form of the coherent matter contribution depends 
in general on the representation of the atomic operators; 
since we are considering identical particles one may use 
the boson representation Ajj. = b^jbk, where [bkjbj] = 
6ki with Ski the Kronecker symbol. The coherent state 
for the matter contribution that preserves the number 
Na of atoms is then written as 


where we have written ajk = Rjk and 7 fc = Qk 
for j, k = 1, 2, ..., n, with gi = 1 and 4>i = 0. Addi¬ 
tionally, we have defined s]- 

An estimation of the ground state energy is obtained 
by minimising the expression (12). Differentiating with 
respect to the phases 9jk and (/)/j one finds critical points 
at 


I 7 ) 





where the Pt operator is defined as 


djk = 0, TT, (j)l - (S'j = 0, TT , (13) 

and by simple inspection the condition for a minimum 
reads 


r'f = 


71 + 72 + ■ ■ ■ + 7n bj 

(|7iP + l72P + --- + |7nP)'/^’ 


( 8 ) 


and where 7 = (71, 72, • ■ ■, 7n)- Here \ 0 )m represents 
the vacuum state in the Fock basis. From the relationship 


[bfc,rt] 


Ik 


A: = 1 , 2 , ..., n, 


one finds the identity [F, Ft] = 1. It is worth noticing 
that this identity shows that F and Pt obey a bosonic 
algebra, hence the coherent state defined in Eq. (7) is 
normalized. One may obtain immediately 




Nalk 


(Er=ii7zP) 


1/2 


[Ft] 


iVa-1 


(9) 


Without loss of generality, one may diminish the num¬ 
ber of variables by choosing 71 = 1. In fact, this assump¬ 
tion is strictly equivalent to eliminate a global phase in 
the matter coherent state (7) and renormalize the other 
values of 7 ^,. 

In order to find the expectation value of the Hamil¬ 
tonian, the matrix elements of the matter operators are 
required. These read 


( 7 l^,fel 7 ) - ^ ^ 

where we have used explicitly in the denominator the fact 
that I 71 P = 1 . 

Finally, the test state that possesses the contribution 
of matter and field is written as 

|a, 7) = |a) 0 I7). (11) 

Using this as a variational state one finds the energy sur¬ 
face 


£ = (a,7|iF|a,7) 


n n 

= £ljk Rjk + A^a Wj 

3<k j = l 



—4-\//W k-jk Rjk 

j<k 


Qj Qk cos{9jk) cosjcSk - (j^j) 

I + Rl 


fijk cos {9^jk) cos {(j>l - <l)j) > 0 . 


Assuming g,jk > 0 the minimum is obtained at 0]k = <t>l- 
(jfj = 0, TT, and replacing it into Eq. (12) the minimum 
energy surface is reduced to 


£ = Y,£l,kR% + N,Y.^j 

j<k j=l 



n 



j<k 


Qj Qk 
1 +^' 


(14) 


Differentiating the energy surface with respect to the 
variables Rjk (related to the expectation value of photon 
number in the jk mode) one finds critical points at 


nc _ r, Q'j Qk _ nrr c 

Rjk - 2 ^kJk 1 + 


(15) 


where stands for the critical value of gk (vide infra). 
Also, since the energy surface is an extensive quantity it 
makes sense to normalize it with respect to the number 
of particles E = S/Na- Additionally, without loss of 
generality we choose wi = 0 , and hence the energy surface 
per particle maybe written as 


E = Y,n,kTfk+Yl 




j<k 


j<k 


i=2 

,c Qj Qk 
l + Rl 


1 + 


(16) 


Finally, the critical points for the variables gk are ob¬ 
tained by solving the system of coupled equations 


dE 

dg, 


Qj 


feU 


1 + ^0 


UJ^ 


-E 


1 + ^0 


-4E 

k=j + l 
1-1 
-4E 


i^jkQk 


Hj/c(l + Rq) 


Q-liQl 


-h 


k^i 
n 

E 


2 22 

Q'j'k'Qj'Qk' 


j'<k 


, D,vfc,(l+i?2)2 


= 0, (17) 


frl 


( 12 ) 
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for j = 2, 3, ..., n, where {p^} = , Q^) repre¬ 

sents the set of critical points. At these points the energy 
surface becomes, using (15), 




r.c 2 


i=2 




n 2 

_4 y ^ 

l<j<k 


i+i?g2 


(18) 


Notice that for any value of the dipolar strengths, g'^ = 
0 and g'j = oo are critical points. In particular, when 
gj = 0 for all j the energy surface takes the value E = 0. 
On the other hand, when a particular g'^ — 0 for s = 
2, 3, ..., n, all the dipolar coupling interactions terms 
related to this level are zero. Thus it is straightforward 
that the atomic part reduces effectively to an {n— l)-level 
system, whose critical points must be analyzed. Hence, 
in this reduction scheme, a given atomic configuration 
will usually lead to different 2-level reduced set; one then 
selects the reduced set with minimal energy in order to 
study the variational ground state of the original system 
If for instance, we consider the limit g 2 —t oo, it is 
straightforward that the energy surface takes the value 


lim E = UJ 2 , (19) 

22 —>00 


which is greater than zero. The interaction terms, which 
give a negative contribution to the energy surface, are 
different from zero (in this limit) if and only if we assume 
that the critical points satisfy 


Qj £>2 ^ 00 , (20) 


for j = 2, 3, • • • ,n. From the expression of the energy 
surface one gets 






i=2 




-4 E 




n 


jk 


Ijlk 




2<j<k 

where 772 = 1 and we have used the relations 

, ^J-ij Vj 


2 


1 + J2k=3 


Q^0'j 




i + EL2^^f i + ELs^r 


( 21 ) 


( 22 ) 


Comparing (21) with (18) it is clear that we have ob¬ 
tained an equivalent system with (n — l)-levels, where the 
atomic variables gj are replaced by the new variables rij 
which upon hnding their critical values lead to an equiv¬ 
alent algebraic system of (n — 2)-levels and so on until 
we reach a two-level system as described above, with one 
radiation mode whose properties, including the complete 
structure of the phase diagram, have been studied exten¬ 
sively [23-26]. 


V 


®3 - T— “3 

^23 ,, M^23 



“1 


FIG. 1. Configurations of a 3-level atom interacting with 
two-modes of electromagnetic field. For the H-configuration 
pi 3 = 0, for the A-configuration pi 2 = 0, and for the V- 
configuration 7733 = 0. The mode of the electromagnetic field 
fljk for each transition is indicated. 


The complete solution of (17) in general requires nu¬ 
merical techniques. However, by using the physical con¬ 
siderations mentioned above, one may find the minimum 
energy surface and hence the variational ground state. 
This new methodology will be exemplified by consider¬ 
ing systems of 3- and 4-level atoms in the next sections 
and the results are generalised for n-level atoms. 


IV. 3-LEVEL ATOMS 


Taking n = 3 in the expression for the variational en¬ 
ergy surface (18), one has 


E = 


1 + Rl 


i^U}2 £?2^ + OJ3 03^ 

M?2 £*2^ 


£^23 £*2^ £*3^ 


{i + Riy 

M?3 03^ 


n 


23 


(l+i?2)2 (l + i?2)2 Oi3 ’ 


(23) 


where without loss of generality we have taken loi = 0. 
In figure 1 the three different atomic configurations are 
shown; also the value fXjk = 0 that dehnes each configu¬ 
ration is given. From the analysis of the previous section 
£* 3 ) = ( 0 > 0 ) is a critical point which yields an energy 
surface equal to zero, E-^ = 0 which defines the normal 
phase. 

By considering = 0 and £>2 7^ 0 j tti® energy surface 
takes the form 


E = 


1 c2 4 

l + gc2‘^2e2 (l + ge2)2 


Hi 2 


(24) 


whose critical points and corresponding minimum can be 
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obtained exactly as 


dividing the parameter space into regions, by 


Q2 = 


I ~ ^12 ^2 

4/i22 + L02 


E-]2 — — 


( — 4 ^12 + ^12 ^ 2 )^ 
16/rf2 ^12 


where the solution exists only for values of the dipolar 
strength 4 /r ^2 ^ ^12 ^ 2 - The equality in this last expres¬ 
sion establishes the border between the normal and a 
superradiant monochromatic region with photonic mode 
^ 12 - 

By considering £>2 = 0 and ^ 0, the energy surface 
takes the form 


E = 


1 


1 + ei 


c 2 

2^3 03 


Ai?3 (?3^ 


(l-te§ 2)2 


(26) 


Its critical point has the same expression as in the pre¬ 
vious case by replacing uj 2 —>■ W 3 , p -12 —>■ /ri 3 , and 
Hi 2 —y that is, 


93 = 


/4^^3 - S2i3 W 3 


4^13 


013 UI3 


Ei3 — — 


("4^13 + ^ 13 ^ 3 )^ 

16/rf3 r2i3 


(27) 

This solution exists only for values of the dipolar strength 
4^33 > f 2 i 3 W 3 , where the equality fixes the border be¬ 
tween the normal and the superradiant monochromatic 
region with photonic mode r 2 i 3 . 

By making the replacement P 3 —> 773 Q 2 into the expres¬ 
sion for the energy surface, and taking the limit g^ ^ 00 , 
one gets 


-E-oo — 


3-^(w2+a;3 773) - 


2 2 
M 23 ^3 


{lE-nlY ^23 


, (28) 


where we identify the energy surface of a two-level system 
with levels given by UJ 2 and 023 . The variable of this en¬ 
ergy surface is 773 . The critical points are 773 = 0, yielding 
an energy value E = 102 , and 


V 3 = 


( 4/^23 ~ ~ ^2)1223 

4/^23 + ~ ^ 2)^23 


(29) 


with corresponding energy given by 


E 23 — UJ2 — 


(4^23 ~ (a^3 — ^ 2 ) ^ 23 ) 
I 6 O 23 M 23 


(30) 


The expressions coincide with the critical energy surface 
of a two-level system when the energy of the lower level 
( 072 ) is different from zero. The latter critical point and 
its corresponding energy are valid in the parameter region 
where the inequality 47 X 23 > ( 0 x 3 — 072)1123 is satisfied. 
Again, the well known expression for the quantum phase 
transition from the normal to the superradiant regions 
for a 2 level system is obtained [23, 24]. 

The phase diagram for the 3-level system is established 
by considering first the equalities 


En = Ei 2 , i?N = Ei 3 , = E 23 , 


4 Mi 2 — ^12 ^2 , 4 7 X 33 — r 2 i 3 073 , 



The first and second equations define bifurcation sets, 
where two different critical points coalesce to (p^j^s) = 
(0,0), while the third equation defines a Maxwell set, 
because the critical points are different but they have 
the same value of the energy. 

We have further divisions of the parameter space, given 

by 


Ei2 — Ei3 , Ei2 — E 23 , Ei3 — E 23 ■ 
The first equality gives 

/ \ 2 

S 2 i 2 


E ?2 = 


X 1 -t 


2 f2i3 7 X 13 
Mi 3 ^13 ^2 


2 ni3 a 73 

Mi3 


± wi 


Mi 3 ^13 W 2 


which must obey the constraints associated to the param¬ 
eters 7 x 12 and 7 x 13 to define a critical point. The second 
equality results in 


2 _ ^23 

7*23 “ 


2 r 2 i 2 /X 32 


2 , ^12 a 72 

7*12 4-^— 


1 + 


7*12 ^12 (*^3 ~ <^2) 


2 if. 


12 


^12 ^2 \2 
4 / 


±1 1 + 


7*12 ^12 (‘X'3 ~ *^2) 

2 j_ r2i2 CJ2 \2 

12 4 / 


(m 


and since we must have 7 X 23 > (073 — 022 )^ 23 /4, only the 
solution with the plus sign is in the collective regime. 
For the third equality one has the first result given above 
by making the replacements 7 x 12 —>■ 7 x 13 , ni 2 —>■ ni 3 , 
u 72 —^ a 73 , and U 73 —^ cj 2 • 

The curves above determine the region boundaries 
where a quantum phase transition between two super¬ 
radiant regions takes place. They correspond to tran¬ 
sitions between the pairs of photonic mode: (f 2 i 2 ,ni 3 ), 
(ni 2 , 023 ) and (ni 3 , 1123 ) respectively. Due to the nature 
of the critical points, they are all Maxwell sets. 

The expectation values for the matter and held ob¬ 
servables with respect to the variational states may be 
obtained at the critical points. In table I we show them 
for the different regions in parameter space. We have la¬ 
belled by Sij the superradian region where the mode Dy 
dominates. In all cases these quantities are proportional 
to the number of particles, and we choose to show them 
per particle. 

By means of the critical points one can give explic¬ 
itly the corresponding variational ground states for the 
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TABLE I. Expectation values and dispersions for the number of photons and atomic populations per particle for the variational 
ground state in 3-level systems, in the different regions of the phase diagram. We have defined (fijk = 4>k — 4’j! s-nd the condition 
for a minimum is (pjk = 0, tt. In the text we have taken oji = 0 without loss of generality. 



different regions in parameter space: 


I^n) = ° ® > 


l^ifc) = 




+ gfcfefc 

. •'/l + , 


Na 


| 0 ) 


M 


(31) 


exp(-iVo r(=fcV2) exp(v^r(=fc |0 )f ,(32) 


Na 


|'I^23) = 


1 


Vsb. 


| 0 ) 


M 


' ^1 + vf 

exp(-W^,j r^l/2) exp( v^r23 “L) lO);’ ,(33) 


where in the second equation fc = 2, 3, and the expres¬ 
sions for r 23 , 02 , £*3 773 are given by equations (15), 

(25), (27) and (29), respectively. The value 7^3 is ob¬ 
tained in the limit 02 00 in the expression for r 23 

assuming 0 § = 03 02 , and yields 


^23 = lim 

Q2^00 


2 fJ-23 V 3 

{1+ r]f )n23 


(34) 


From these general results for 3-level systems we can 
easily obtain the phase diagrams for the different atomic 
configurations. For the S-atomic configuration (7113 = 0) 
one has the critical energy surfaces E^, £’ 12 , £23, from 
which the phase diagram can be obtained and is displayed 
in figure 2 (a). 

For the F-configuration {fJ ,23 = 0) we have critical en¬ 
ergy surfaces £ 12 , £ 13 , and the corresponding phase 
diagram is plotted in figure 2 (b). 
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Finally, for the A-configuration {^12 = 0 ) the mini¬ 
mum energy surfaces are Em, E13, E23, and the quantum 
phase diagram is shown in figure 2 (c). 


In general we may write 



Esmin = minjEN, ^12, E23} , 


m 

Evmin = min{EN, Ei2, E13} , 



E\min = minjEN, E13, E23} , 

( 35 ) 



where Em exists for all values of fiij, E12 is independent 
of /ii3 and 1123, and so on. Other critical values for g2 and 
£>3 exist; however, in all cases they provide higher energy 
values. Notice that since the solutions for E12, E13, and 
E23 only have emission or absorption of photons associ¬ 
ated to the modes O12, O13 and O23, respectively, the 
polychromatic phase diagram for the variational ground 
state is naturally divided into monochromatic subregions, 
labelled by Sij in Fig. 2 . 

In figure 2 (a), the minimum energy surface as a func¬ 
tion of the dipolar intensities /ii2, /i23 is presented. The 
normal region Em = 0 is shown in black, while the 
collective regions, with values E-^min < 0, appear in a 
color graded scale. The separatrix (white lines) is de¬ 
fined by the set of points {gi2,g-23) where E^ = E12 
and FIn = E23, indicating the borders between the nor¬ 
mal and collective regions, together with the points where 
Ei2 = E23 dividing the collective regimen into monochro¬ 
matic regions, corresponding to the two electromagnetic 
modes. 

According to the Ehrenfest classification [ 22 ], a phase 
transition is of order j, if j is the lowest non-negative 
integer for which 


(b) 



(c) 


1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 


1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 


1.2 



d^E^ 


lim „ . 
£->■0 




7S lim ——r- 

e-i-o dfj,^ 


where fj, represents a dipolar intensity parameter. In our 
case the ground state undergoes a second order phase 
transitions when the dipolar intensities go from region N 
to region S'12 and vice versa, but first order phase tran¬ 
sitions occur when crossing from region N to region S23 
and from region 5'i2 to region 5 ' 23 . 

As an example of the behavior of the expectation val¬ 
ues given in Table I, figure 3 (a) shows {1^12) = (a|2ai2) 
and {A]^k) (k = 1, 2) as functions of fii2 for ^,23 = 0; fig¬ 
ure 3 (b) shows (1/23) = (a23®23) and {Akk) {k = 1 , 2 , 3 ) 
as functions of ^,23 for g,i2 = 0 ; and figure 3 (c) shows 
(Vjk) = {a]k^jk) and (Akk) {j < k = 1, 2 ) as func¬ 
tions of (^, where we have parametrised the dipole inter¬ 
actions in the form gi2 = /rcosC, g-23 = /xsin^, /i = 1, 
and 0 < C ^ '^/2, for a trajectory going through re¬ 
gions S'12 and S23. The breakpoints and discontinuities 
inherited from Esmin and its derivatives at the loci of 
a phase transition are evident. For first order transi¬ 
tions the expectation values exhibit a discontinuity. We 
have taken fli2 = 1 , ft23 = 0 . 5 , uj2 = 1 , UJ3 = 1 . 3 . 
The same results are obtained for all values of g,i2 that 
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0.8 

m 
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FIG. 2. (color online) Energies and phase diagrams for atoms 
in the H-, V-, and A-conhgurations. The separatrices (white 
lines) and the order of the transitions are shown. The normal 
regions are labelled by N (black). The collective regime is 
divided into the regions: (a) S '12 and S 23 corresponding to 
the active modes fli 2 and Q 23 respectively, shown here for 
fli 2 = 1, 1123 = 0.5, 012 = 1 and 013 = 1.3; (b) S 12 and S 13 
corresponding to the active modes Q .12 and flis respectively, 
with fli 3 = 1 , Q ,\2 = 0 . 8 , UJ 2 = 0.8 and = 1 ; and (c) 
S 13 and S 23 corresponding to the active modes lli 3 and 1123 
respectively, with II 13 = 1 , II 23 = 0 . 8 , U 12 = 0.2 and CJ 3 = 1 . 
































FIG. 3. (color online) behavior of the expectation values of 
the photon and atomic level populations, as functions of the 
interaction parameters. See text for details. 


satisfy 4/rj2 — 0J2^i < 0, or values of fi23 satisfying 
- (^3 - W 2 ) rj 23 < 0 (see Fig. 2 (a)). 

Figures 2 (b) and 2 (c) show maps of the minimum 
energy surface for atoms in the V- and A-configurations. 
The normal region is shown in black, and the collective 
region in a color graded scale. For the V case we have 
second order transitions when the system goes from N to 
S'12 and from N to S13, and first order transitions when 
it goes from S'12 to S13. A similar analysis of the physical 
quantities (number of photons and atomic populations) 
to that of figure 3 will show that they all vary contin¬ 
uously as functions of fii2 with d/rfg — 0^3 < 0 and 

as functions of p,i3 with 4/if2 — W2 ^12 < 0, because of 
the second order of the transition, while they exhibit dis¬ 
continuities when the system goes from region S12 to S13 
as this is a first order transition. For the A case, we 
have UJ2 > wi, the phase diagram for atoms in this A- 
configuration resembles the case of the S-configuration 
(compare figures 2 (a) and 2 (c)). The expectation values 
of physical quantities vary continuously when the system 
goes from N to 513, where a second order transition oc¬ 
curs, and they present discontinuities when the systems 
goes from N to S23 and from 5'i3 to S'23, where the system 
undergoes a first order phase transition. 

All the results for 3 -level atoms show that the poly¬ 
chromatic minimum energy surface divides itself into 
monochromatic regions where only one mode of elec¬ 
tromagnetic field contributes to the variational ground 


state. Also, as first order transitions are related to dis¬ 
continuities in at least one physical quantity, these criti¬ 
cal points along the separatrix form a Maxwell set, ie., 
the corresponding critical points in the transition = 0 
and —>■ 00 provide the same minimum energy value, 
whereas the critical points where second order transitions 
occur correspond to bifurcations. 


V. 4-LEVEL ATOMS 


Taking n = 4 in the expression of the variational en¬ 
ergy surface ( 18 ) one has 

1 _|_ (^*^2 Q2^ + W3 £>3^ + a;4£l4^^ 


E = 


1^23 £>2^ 03^ 


M?2 £*2^ 


{I + RIY 

4 


^^23 

1 A 3 ef 


(l+i?2)2 O 12 

4 


{l + RlY f2i3 

4 


(l + i?2)2 fl34 


£^24 £>2 ^04^ 


{I + RIY 

4 


O 24 

,,2 „ c 2 c 2 

M34 Qz £'4 


{l + RlY 1234 ■ 

Again, {g^, 93, Q4) = ( 0 , 0 , 0 ) is a critical point which 
yields an energy surface = 0 . 

By considering el = 0 and Q% ^ 0 with k = 2 , 3 , the 
energy surface takes the form of a 3 -level system 


E = 


1 


1 + 4?o 


^a;2 02^ + 1^3 £*3^ 


M 23 


(l + i?2)2 f|23 


M12 £*2^ 


M?3 Bz^ 


{IERif n,2 {i + Rlf ^iz ■ 

By further taking el = 0 and ei = 0 we have 

(—4 £132 + f 2 i 2 ^2)^ 


( 37 ) 


£>2 = 


I 4/if2 ~ ^12 ^2 

4/^12 4” ^12 ^2 


E-12 — - 


while for el = 0 and ei = 0 one gets 

£>3 = 


1 4^13 ~ ^13 ^3 


4/^13 + ^13 W 3 


Eiz — — 


16/4^2 ^12 

( 38 ) 

(— 4/433 + ^13 ‘^zf 


16/4^3 Oi3 


( 39 ) 

where the previous expressions exist only for values of 
the dipolar strength 4 /t^^ > ftik ujk, with k = 2 , 3 . 

In the limit when the variables go to infinity, one has 
el = 0 , ei oOj Bz — Vz £’21 the critical point given 
by 


Vz = 


' 4/423 ~ ('^3 ~ <^2)023 


4/423 + ('^3 ~ <^2)1223 
ly surface takes the ex; 
(4/423 ~ (*^3 ~ ^ 2 ) ^2z) 


( 40 ) 


from which the energy surface takes the expression 

E2Z = UJ2 — - -7777;;-2- > ( 41 ) 


161123 £423 
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^^34> ^34 


[l23, il23 


M^12> ^2 



<A44> = 0 


^23’ ^23 


h2’ ^12 




FIG. 4. (color online) Under the condition (A 44 ) = 0 for the 
ground state, a system of 4-level atoms in the H-configuration 
is reduced to 3-level atoms in the H-configuration. 


which is valid in the parameter region where the inequal¬ 
ity 4/i23 — (^3 ~ ^ 2 ) 1^23 is satisfied. The expression 
coincides with the critical energy surface of a two-level 
system. One again, the well known expression for the 
quantum phase transition from the normal to the super¬ 
radiant regions for a 2 -level system is obtained. 

4-level systems have many configurations named by 
the atomic physics community by H, A, y, N, <(>, ..., and 
the number of modes yielding transitions between pairs 
of atomic levels depend on the considered configuration. 
For 4-level atoms in the H-configuration, to take q1 = 0 
is equivalent to decouple the fourth level from the sys¬ 
tem and the population of this level therefore vanishes. 
Thus the system is reduced to 3-level atoms in the same 
configuration as shown schematically in hgure 4. The 
A- and ^-configurations, under the same condition, are 
reduced to the 3-level atomic configurations A and V 
respectively, as shown in figure 5. On the other hand, 
under the condition P 2 = 0 (to be considered below) the 
A- and <)>-configurations are reduced to 3-level atomic H- 
configurations as shown in the same figure. 

By considering the case = 0 and g% ^ 0 with k = 
2, 4, the corresponding energy surface takes the form of 
a 3-level system by making the replacement of the label 
3 —)■ 4 in equation (37) 



FIG. 5. (color online) Schematic procedure to obtain, by a 
simple consideration {{Akk} = 0), a 3-level systems from a 
4-level system. The A-configuration is reduced to the H- or 
A-configuration and the <)>-confignration is reduced to the H- 
or U-configuration. In both cases the simple condition on the 
ground state is given by {A 22 ) = 0 and (A 44 ) = 0 , respec¬ 
tively. 


for values satisfying 4^44 — fli 4 a ;4 > 0 . 
In the limit ^ 00 and g1 = 774 P 2 j 


Vi = 


' 4^24 ~ (^4 - 072)024 
4/724 + - ^2)^2i ’ 


(45) 


W 3 —t 074 , r2i3 —> f2i4 , 1223 O 24 , 

03 P 4 , /7i3 —)■ /7i 4 , /723 —t g,2i ■ (42) 

For the case 03 = 0 and 04 = 0 one finds the result 
already considered above. Then if we take = 0 and 
P 2 = 0 , the critical points and energy surface are given 
by 


/ 4/7^4 - 1214 074 
y 4/7 i4 -I- 12i4 074 

(43) 

( —4/7i4 -|-12i4 074)^ 

-^14 1 2 0 ’ 

i2i4 

(44) 


and the energy surface is 

(4^24 ~ (^4 - OJ 2 ) ^ 24 )" 


E2i — UI2 — ' 


16 O 24 uli 


(46) 


with 4 / 7 I 4 — {uji — 072)1224 > 0. 

Finally, for g^ = 0 and 0 ^ 7 ^ 0 with k = 3,4, the 
corresponding energy surface takes the form of a 3-level 
system with the replacements of the labels 2—7 3 and 
3—7 4. Then if we consider g\ = 0, or pg = 0, we have 
the results already discussed above, i.e., E 13 and En. 

In the limit when the variables go to infinity, one has 
02 = Oj 03 00 , gl = fji g'^, and thus the critical points 
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are given by 


Vt = 


l ^fJ-34: - (^4 - ^ 3)^34 

4/i§4 + (a;4 - a;3)034 


(47) 


with 


Esi = ^^3 — 


(4^34 - i.^A- ^3)^3i)^ 

16 O 34 


(48) 


when 4 — (a ;4 — a; 3 )ll 34 > 0 is fulfilled. 

The phase diagram for a 4-level system is established 
by considering the equality of the different energy sur¬ 
faces, that is, from 



TpC _ TT'C TTiC _ TT’C TT'C _ 77’C 

£/N — -^12 5 — -^13 7 -^N — -^23 ’ 

El^ = El^, i??2 = i?23, ^^?3 = ^^23, (49) 

giving the expressions obtained already in the 3-level 
case. Therefore we will have a separatrix formed with 
bifurcation and Maxwell sets. 

The remaining conditions are associated to the new 
critical energy surfaces by establishing the equalities 


FIG. 6. (color online) Phase diagram as a function of the 
dipolar intensities for atoms of 4-levels in the H-configuration. 
Parameters are: fli 2 = 1, ^23 = 0.7, 1134 = 0.3, cji = 0, 0^2 = 
1, 013 = 1.7 and u }4 = 2. 


VI. n-LEVEL ATOMS 


£& = 

= £f4, 

£n 

II 

£n = 

0 ^ 

II 

£?2 = 

= £f4, 

EI 2 

= £24 , 

EI 2 

= £ 3 ^ 4 , 


= £m : 

EI 3 

= E^24 : 

EI 3 

= £ 3 %, 

£23 = 

= £m , 

£23 

= E^24 , 

£23 

= £ 3 %, 

n 

= E^24 , 

£f4 

II 

£24 

= £34 • 


The previous results determine the total number of bor¬ 
ders associated to 4-level atoms, however some of them 
are not allowed by the selection rules of the dipolar in¬ 
teraction. To see this, we consider the S-configuration 
where one has /ii 3 = /ii 4 = ^ 2 A = 0. Then the quantum 
phase diagram is divided into 6 regions separating the 
normal and collective regimes 


We may generalise the above results to n-level atoms 
in dipolar interaction with (. modes of electromagnetic 
field. 

We have found three types of minimum energy sur¬ 
faces: 


• Epf = 0, corresponding to critical points Pg = 0 for 
s = 2, 3,..., n. 


• Elk, when flik promotes the transition wi ^ uik 
with Pg = 0 for all s ^ /c, given by 


^ (Wfc flu. - 4/^2^)^ 

^Ik — A a 2 

16 


(52) 


TT'C _ TT'C 771C _ TT’C TT’C _ TT'C 

£/N — -^12 5 — -^23 7 -^N — ^34 5 

^12 = E23 , E-^^2 = -E'34 , ^^23 “ ^34 : (^9) 

whose explicit expressions and existence condition were 
established previously. 

The minimum energy surface for the this configuration 
of 4 levels is given by 

E’^min = minjilN, £^12, £23 , £34} , ( 51 ) 

where, again, one must take into account the regions in 
parameter space where they are valid. In figure 6 the 
separatrix of this system is shown. Clearly the energy 
surface is divided in a normal region and three collective 
regions, where in each of the latter only a monochro¬ 
matic electromagnetic field contributes to the ground 
state. One finds second order transitions when the sys¬ 
tem goes from region N (normal region) to region S' 12 , 
while for other crossings the transitions are of first order. 


when the condition 4 — ujk ^ik > 0 is fulfilled, 

with the remaining critical point being 


Qk = 


UJk fllfc 


(53) 


• Ejk, when fljk promotes the transition ujj ^ uJk, 
with = rjc in the limit —>■ 00, given by 


Ejk — 


(^{uJk-ujj)njk-Aji^^^'^ 


16 0 




(54) 


when the condition 4/4^^ — {ojk — ojj) fljk > 0 is 
fulfilled, with the remaining critical point rj^ being 


Vk = 


I'^t^jk- i^k - kkJj) ^jk 

Ag?-k + (wfc -u}j)njk’ 


(55) 
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From the above, the minimum energy that of the vari¬ 
ational ground state is given by 

Emin = minlFl^v, {Eik}, {Ejk}} , (56) 

for j < fc = 2, 3, ..., n. 

A. Separatrix and order of transitions 


B. Expectation values 

The expectation values of matter and field observables 
for the variational ground state, such as number of pho¬ 
tons and atomic population, may be calculated analyti¬ 
cally. That for the product of matter and field observ¬ 
ables reduces to the product of the expectation values of 
each one, i.e., 


The separatrix divides the minimum energy surface in 
regions where the properties of the ground state change 
significantly. An example of this is the boundary be¬ 
tween the normal region, where the ground state is given 
by the vacuum state for the field and all atoms in their 
lowest level, and the collective region where the ground 
state possesses contributions from states with a nonzero 
number of photons and/or excited atomic states. 

For values of the dipolar intensities satisfying 


^ ^ik , 


4 , 


(57) 


the minimum energy value is Emin = Ejy = 0. The 
separatrix that divides the normal and collective regions 
is given by the points where the equality in equations (57) 
is attained. The separatrices in the collective region are 
given by the sets of equalities Ejf^ = Eji^i, for the modes 

^jk ^ ^j'k'- 




where 


G = + i^k' - 

8 ^j'k' (2 UJi UJk' ) 5 

with i = j, fc, and the appropriate sign is taken to lie in 
the collective regime. 

From the expressions for the minimum energy value in 
the different regions, and following the Ehrenfest classi¬ 
fication [22], one finds the order of transitions at each 
separatrix as 


Transition 

Order 

En ^ Elk 

Second 

En ^ Ejk 

First 

^jk Ej'k' 

First 


We stress that the critical points where a second order 
transition occurs form bifurcations, i.e., where the points 
coalesce at the origin in the p-space and hence all physical 
quantities (expectation value of the number of photons, 
atomic populations, etc.) change in a continuous manner. 
The critical points where a first order transition occurs 
form a Maxwell set, i.e., there are two different critical 
points with the same energy value. In this case at least 
one physical quantity has a discontinuity. 


(Om <8) Of) = {Om) {Of) ■ (60) 

For this reason, in this variational approximation, there 
are no correlations between matter and field operators. 
In a similar manner as with the minimum energy surface, 
it is necessary to consider the normal and the different 
collective regions. In the normal regime one finds that 
the variational ground state is given by 

IV'w) = ® 1^)^’’ (®1) 

where |0)m and |0) f denote the vacuum state in the Fock 
basis for matter and field respectively. Clearly, we have a 
zero expectation value for the number of photons, and the 
atomic populations for this state is the corresponding to 
the lowest atomic level, i.e., (An) = 1 (per particle). On 
the other hand, in the collective region, if the minimum 
energy surface is given by Eji^ for fixed values of j and 
k with j < fc = 2, ..., n, the variational ground state 
takes the form 



0 exp(-7VarJ^/2) exp(\//W aj^,) |0 )f , (62) 
when J = 1, with 

Q% 


Tlk = 


(1 + Q‘i^)^lk 


(63) 


and 




1 lb] + vlbl 


No. 






| 0 ) 


M 


0 exp(-7Va f/l/2) exp(v^f/fe ajj |0 )f , (64) 
when J > 1 and with 


r 


c _ 

jk — 


lim 


'jk 


e%=v% Bj 


(1 + Vk'^)^jk 


(65) 


The indices jk indicate that we only have photons of 
mode Qjk and that only the atomic levels j and k are 
populated, with zero contribution from the rest. Ta¬ 
ble II, shows the different expectation values and fluc¬ 
tuations for the normal and collective regions Sjk- One 
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TABLE II. Expectation values and dispersions of the number of photons and atomic populations per particle for the variational 
ground state in n-level systems, in the normal N and collective regions Sjk of the phase diagram. We have defined (pjk = cj>k — 4'jj 
and the condition to have a minimum is cpjk = 0, tt. In the text we have taken oji = 0 without loss of generality. 



should note that the photon distribution is a Poisson dis¬ 
tribution, i.e., one has (A UjkY' = (I'jk), while the matter 
contribution is described by a binomial distribution 

(^“) (66) 

where p stands for the atomic population per particle of 
the jth level and q for the corresponding one of the A:th 
level, i.e., from table II one has 

_ 1 / (gjfc ~ \ 

H “"k J ’ 

and hence the corresponding atomic populations are re¬ 
lated to them by {Ajj) = NaP, and (Akk) = Na q. 


C. Symmetries in the RWA 

When the RWA approximation is considered the terms 
that do not preserve the total number of excitations in 
the Hamiltonian are neglected, i.e., the Hamiltonian is 


written as 


= Hrt — 




n 

'y ( Mjfc {^\k A -^kj^ ) (67) 


j<k 


where Hd is given by Eq. (2). Using the same test state 
discussed previously, one finds the corresponding varia¬ 
tional energy surface in RWA approximation to be 


^RWA = ^ R^jk + A^a ^ W 
j<k j=l 


^ 1 + e :=2 9 i 


y ( P’jk Rjk 

j<k 


Qj Qk COs(()>fc 

1 +ru 0i 


Comparing with Eq. (12) one sees that the effect of 
considering this approximation modihes the dipolar in¬ 
tensities by a factor of 1/2; we may therefore find the 
solutions to the problem in this approximation by re¬ 
placing pLjk —>■ Pjk/‘i in all expressions for the critical 
points of the complete Hamiltonian. Concerning the 
phases, the critical points are given by the relationship 
'('fc ~ 4’j ~ ^jk = 0, TT with the condition for being a min¬ 
imum fXjk cos{(j)l - (j^j - (^jk) > 6- 

On the other hand, one may prove that the Hamilto¬ 
nian (67) possesses n linearly independent constants of 
motion Kj (including the total number of particles Na), 
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TABLE III. 3-level atomic configurations: shown are the dif¬ 
ferent values of Xk which correspond to the number of excita¬ 
tions that are required to excite a single atom from its lowest 
level to fcth level. 


configuration 

A 2 

A 3 

s 

1 

2 

V 

1 

1 

A 

0 

1 


by using the relationships 

ag{i^) = g{v+ l)a, ( 68 ) 

gii') = g{i^ - ^) a.\ (69) 

■^jk gi.-^jji — g^-^jj ^kk T 1 ) X^jk ; (^ 6 ) 

^jk g^^jjt -^kk') ~ g(.-^jj ^5 X\.kk 1 ) ^kj ; (^ 1 ) 

for field u = a and matter Ajk = b^- operators, and 
where g{-) stands for an arbitrary analytical function. 
These constants of motion are given in general by 

= -^n + ^kj — i^jk ■ (72) 

k<j j<k 

Clearly the first order Casimir operator in terms of these 
constants is given by 

n 

Na=Y.K,. (73) 

Also, the total number of excitations is written as 

n n n 

M = ^ Xg Ki = ^ Ujk + ^ Xk Akk (74) 

t=2 j<k k=2 

where the integer values of Afc, which depend on the par¬ 
ticular atomic configuration, stand for the number of ex¬ 
citations that are required to excite one atom from its 
lowest atomic level to the /cth atomic level. As an exam¬ 
ple of this table III shows the values of Xk for the case of 
3-level atoms in their three atomic configurations. 

It is interesting to note that for each constant of mo¬ 
tion of the Hamiltonian in the RWA approximation (67) 
one has a corresponding symmetry operator for the full 
Hamiltonian (1). This is given by 

Ilj = exp {i TT Kj). (75) 

These operators may be useful to define the symmetry- 
adapted variational states, which provide a better ap¬ 
proximation to the exact quantum ground state. 


The constants of motion in the RWA approximation 
will lead to a better variational test function as was 
shown for the two-level atoms [24, 27]. This is done 
by truncating the power series expansion of the coherent 
states. For the generalized Dicke case the parity sym¬ 
metry operators will be useful to define the symmetry- 
adapted variational states, which provide a better ap¬ 
proximation to the exact quantum ground state. 

VII. COMPARISON WITH THE QUANTUM 
GROUND STATE 

In order to compare the above variational solution with 
the exact quantum one, we consider 3-level atoms in the 
S-configuration. Using the Fock basis 

Wi2, 1^23, ni, n2, ns), (76) 

for the two modes and the three atomic levels one may 
diagonalise numerically the Hamiltonian in order to find 
the ground state. However, one may use the fact that 
this system possesses three different symmetries Ilj = 
exp^in Kj) (75) for j = 1, 2, 3 with 

Ki = All — 1^12, 

K 2 = A 22 + *>'12 — *^ 23 , (77) 

Ks = A33 -I-1/23. 

Also, both the total number of excitations and 

the number of atoms Na may be written as a linear 
combination of them as ATh = K 2 + 2 Ks and Na = 
Ki -I- K 2 X- Ks- Using these and Ks one may rewrite the 
Fock basis as 

\M-Na-Ks + ni, Ks-ns, ni, Na-ni-ns, ns ), (78) 

where M and Ks are the eigenvalues of A4 h and Ks: 
respectively. In fact, using the corresponding symme¬ 
try operator for the total number of excitations 11 ^ = 
exp{iTrMs) and the 1 X 3 operator one finds that the 
Hamiltonian (1) is divided into four blocks which pre¬ 
serve the parity of the and Ks operators, i.e., the 
Hilbert space is divided into four ortogonal subspaces 

{IV'ee)}, {\lpeo)}, {IV'oe)}, {|'0oo)}, (79) 

where the subindexes e = even and o = odd denote the 
parity of ATh and Ks, respectively. Using these subbasis 
the Hamiltonian may be diagonalised numerically and 
hence the minimum energy value E^e, E^o, Egg, E^o for 
each subspace may be found. Thus, the exact quantum 
ground energy is given by 

Emin = min{A ee; Eeo, E oe 1 Eoo} ■ (80) 

Figure 7 shows the quantum ground energy per par¬ 
ticle (80) (color graded scale) in comparison with the 
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FIG. 7. The quantum ground energy surface Emin (color 
graded scale) in comparison with the corresponding varia¬ 
tional ground energy (mesh surface) as a function of the con¬ 
trol parameters for a number of particles (a) Na = 1 and (b) 
Na = 2. 


FIG. 8. The rate of the difference to the sum of the expecta¬ 
tion value of the number of photons, Sv, of the exact ground 
state for (a) Na = 1 and (b) Na = 2. 


corresponding variational ground energy (mesh surface). 
One observes in both cases Na = 1 Fig. 7(a) and Na = 2 
Fig. 7(b) that the variational ground energy surface is a 
very good approximation to the exact quantum calcula¬ 
tion. 

On the other hand, it is of interest to investigate the 
polychromatic behavior of the phase diagram. For this, 
we consider the ratio of the difference to the sum of the 
expectation values of the number of photons, defined by 


{ i >23 - t^l2) 
{U 23 + 1^12) 


(81) 


which takes values —\<5v<l except when the state 
satisfies {h >23 + i>'i 2 ) = 0 , which is satisfied only in the 
normal region with ^12 = 0 (black region in figure 7). 
Clearly, when <5 « — 1 the ground state is dominated by 

the mode O 12 , in contrast with values 5v where the 
state is dominated by the mode fl 23 . 

Figure 8 shows the ratio of the difference to the sum of 
the expectation values of the number of photons, (5 v, for 
the corresponding ground states of figure 7. The classical 


separatrix is shown in white lines. One may note that this 
quantity is not sensitive where second order transitions 
occur (cf. Fig. 2 (a)) and takes values 61 / Ki —1 in both 
the normal and collective regions, indicating that the 
mode O 12 dominates, but around the separatrix where 
the first order transition occurs both modes f 2 i 2 and O 23 
contribute to the ground state, since Sv^O. The region 
above the separatrix is dominated by the mode O 23 where 
one has Jz/ « 1. We note that the quantum separatrix 
where a first order transition occurs {Si' = 0 , black re¬ 
gion in figure 8 ) approaches to the variational separatrix 
when the number of particles grows, in fact for Na = 2 
the quantum separatrix is very close to the variational 
one Fig. 8 (b), in contrast with the case Na = 1 Fig. 8 (a). 

The expectation value ( 2 ^ 12 ) is negligible in the normal 
region; this fact is shown in figure 9, and we see that the 
quantum phase diagram is divided into monochromatic 
regions in a similar fashion to the variational calculation. 
This is now only in the sense that the bulk of the ground 
state is dominated by states with only one type of pho¬ 
tons in each region; in fact this figure shows that in the 
normal region the ground state is dominated by the vac- 
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FIG. 9. Expectation value of the number of photons lyw) of 
the exact ground state for (a) Na = 1 and (b) Na = 2. The 
classical separatrix is indicated (white line). 


uum state. 


VIII. CONCLUSIONS 


A system of Na atoms of n-levels interacting dipolarly 
with £ modes of electromagnetic field, where the transi¬ 
tions between two given atomic levels are promoted only 
by one mode, has been considered. Using as a variational 
test state the direct product of Heisenberg-Weyl states 
for the field contribution, and coherent states preserv¬ 
ing the number of atoms for the matter contribution, the 
variational energy surface is calculated and the minimum 
obtained. 

We have given an iterative procedure for the calcu¬ 
lation of the critical points that reduces any system of 
n-levels to 2-level systems. 

The procedure to find the critical points that min¬ 
imise the energy surface was exemplified for 3- and 4- 
level atoms interacting with ^ modes {£ = 2 for 3-level 
atoms and £ = 3, 4 for 4-level atoms). The normal and 
collective regions were described analytically and was 


demonstrated that the collective region is divided into 
I monochromatic regions where only one mode of elec¬ 
tromagnetic field contributes to the ground state while 
the other ones remain in the vacuum state. 

Studying the transitions between the different regions 
in the phase diagram, for 3-level atoms, we find both first 
and second order quantum phase transitions (cf. Fig. 2). 
First order transitions are directly related to the fact that 
at least one physical quantity of the system changes in 
a discontinuous manner (cf. Figs. 3(b) and 3(c)), which 
is related to the fact that the set of critical points in the 
separatrix forms a Maxwell set. Second order transitions 
present a continuous behavior (cf. Fig. 3(a)), and in this 
case the critical points form bifurcations. Similar results 
are obtained for the general case of n-level atoms inter¬ 
acting with £ modes. 

This variational study suggests the following relation¬ 
ship between matter and field observables for the ground 
state: 

( 82 ) 

which, at the separatrix and considering wi = 0, implies 

(Ayt„)V (83) 

i Ljf. 

Being this a universal relationship, one may propose it 
as an experimental criterion to detect the transition be¬ 
tween the normal and superradiant regimes. 

We have shown that, when the RWA approximation 
is considered, the phase diagram of the system suffers 
a rescaling of the dipolar intensities by replacing /r —>■ 
(/j,)'^^'^/2 in all quantities. Also, we have shown that 
in the RWA approximation the Hamiltonian possesses n 
linearly independent constants of motion (including the 
total number of excitations), and that each one of them 
provides a symmetry operator of the full Hamiltonian. 
These will be useful to establish symmetry-adapted vari¬ 
ational states. 

Finally, using the symmetries of the full Hamiltonian 
we found by numerical calculation the exact quantum so¬ 
lution for the case of 3-level atoms in the S-configuration. 
We also found that the variational minimum energy sur¬ 
face is a very good approximation to the exact quantum 
one. Also, we considered the ratio of the difference to 
the sum of the number of photons and found that the 
quantum phase diagram is divided into monochromatic 
regions, this in the sense that the bulk of the ground 
state contains contributions of only one type of photons 
in each region. 
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